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ABSTRACT 

This aper reviews the information relating 
the average* sustainable yield (Y) from a fishery 
resource to the bicnass (BQ) present before ex- 
ploitation starts. By examining the population 
parameters of a range of fish stocks, and using 
simulation modelling it has been shewn that the 
cotitonly used formula Y = 0.5 M BO gives too 
high a value for Y. Hie degree of over-estimation 
depends on the nature of the stock-recruitment 
relation, and on the variability of recruitment. 
The iitplications for nanagenent and survey design 
are discussed. 
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1. INTRODUCTION 

The background of this study is concerned with the problem of estimating the poten- 
tial sustainable yields of fish stocks that hitherto have been unharvested. The need 
for guidelines for such estimation of yields is obvious. However, this need has been 
exacerbated by the extension of exclusive economic zones and the considerable interest 
that this has engendered in coastal states in calculating the potential of their newly 
acquired resources. 

In responding to these demands FAO has been involved in a series of surveys using 
the R/V Fridtjof Nansen (FAO 1978; Ke Steven, Nakken and Stromme 1981) and other vessels 
(Troadec and Garcia 1980; FAO /UNDP 1981; Vidal-Jtlnemann 1981). These surveys use a 
mixture of acoustic and other techniques to produce estimates of the abundance of a 
particular fish stock. An approximation to the maximum sustainable yield (MSY) of this 
biomass is then calculated using some simple formula. The most commonly used form is 

MSY 1/2 MB 
o 

where M is the coefficient of natural mortality and B the estimated biomass of the fish 
stock. A variety of other forms have been proposed to take account of previous levels 
of exploitation and yield curves which reach their MSY at other than 1/2 B (Gulland 1971). 
This approximation to MSY is then used as an indicator of the potential yiild of the 
resource. 

In this study we examine the validity of this approximation both in general, and in 
particular, where it has been used when coupled with research surveys. 

It is a beguiling idea that if one knows sufficient of the life history character- 
istics of a species, it will be possible to make some estimate of its potential for sus- 
taining catches. It is an even more beguiling idea that there is some simple expression 
which encapsulates that information and permits sustainable yields to be calculated for 
all species. This in essence is the idea of the so-called production model which sub- 
sumes all the details of life history characteristics and density-dependent responses 
into the single equation in biomass B. 

For the Schaefer model the particular form used is 

-'B(1-B/B o > (1) 

Here the MSY level occurs at half the virgin biomass, B , and this feature, coupled with 
the claim that maximum sustainable yields appear to occur when fishing mortality F is 
similar to natural mortality M, leads to the use of the expression 

MSY - 1/2MB (2) 

o 

as an approximation to the maximum sustainable yield. 

Recently Shepherd (1982) has made some similar calculations for more general produc- 
tion models and concluded that the MSY tended to be somewhat less than the expression 
1/2 MB , but depending on choice of parameters greater values could be obtained. 

This use of the expression presupposed that the production model is a reasonable 
representation of the dynamics of fish stocks. Some check on this can be made by con- 
sidering equation systems of a more realistic nature. For example a more detailed 
model may be constructed using the framework of the Beverton and Holt equations and 
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estimates of sustainable yield can be made for different parameter combinations. Gull and 
(1971) adopted this approach and reached the broad conclusion that within a reasonable 
range of parameter values and fishing strategies the use of equation (2) would produce a 
satisfactory approximation to the maximum sustainable yield predicted by the "simple 11 
(constant recruitment) Beverton and Holt model. This result has led to the use of ex- 
pression (2) in calculations of the potential yield of a wide variety of fish stocks 
(Gulland 1971). 

There are three major problems with the use of this approximation. Firstly, the 
analysis of the Beverton and Holt model, considered above, is misleading in that for most 
realistic choices of an age or size at recruitment the MSY as a proportion of unexploited 
biomass is rather less than the 1/2 M level. This problem of choice of age at recruit- 
ment had been mentioned by Gulland (1971), but was not fully explored by him. 

Secondly, harvesting at the MSY level in some cases involves a reduction of the 
spawning stock biomass to a level where the chance of recruitment declines are signifi- 
cant, hence the assumption of constant recruitment would be likely to be violated. One 
possibility, explored below, is to constrain the harvest level to ensure that there is 
some constant escapement. 

Thirdly, as is well known, recruitment is highly variable in some fish stocks even 
when the parent stock is relatively unchanging. This last problem has several impli- 
cations. Firstly, even though harvesting at MSY may be in some cases be compatible with 
a given escapement if recruitment were constant, recruitment variation may mean that 
there is a high probability that there will be a reduction in spawning stock to below 
this escapement level. Secondly, a quota management system which seeks to harvest a 
fixed quantity each year, corresponding to the estimated MSY, may lead to unacceptably 
high variations in fishing mortality and effort. Thirdly, a target effort level based 
on MSY may produce unacceptably high variation in catch. Fourthly, survey estimates 
of initial biomass may be misleading as the biomass fluctuates with recruitment, and 
the surveys may have been made when abundance was uncharacteristically high or low. 

The remaining analysis is organised into four sections. In the first the basic 
Beverton and Holt model is re-examined, and the actual MSY levels for different para- 
meter combinations are presented. 

In the second, the effect that MSY harvesting has on the spawning stock is examined 
and the problem of an appropriate escapement level is discussed. 

In the third section the underlying variability of recruitment is considered and the 
fate of model populations that undergo the process of survey, assessment and harvesting 
are explored by Monte Carlo techniques. In this section the potential yield calculations 
are modified to ensure that the probability that a stock will be reduced below some fixed 
escapement in a particular time span is less than some specified value. 

In the final section the results of the analysis are used to suggest a series of cal- 
culations that need to be made to assess the potential of an unexploited fish resource. 

2. THE BEVERTON AND HOLT MODEL REVISITED 

In this section the analysis is developed in a simple form of the Beverton and Holt 
model. 

^ 

X. is defined as the numbers of age class i in year t, w. as the mean weight of age 
class i. R is the recruitment, w and K the parameters of tfte Von Bertalanffy curve, 
and F and M the coefficients of fishing and natural mortality. 
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The dynamics in the unexploited state are then given by the equation system: 



X o,t - R 

X, . = X, - . - exp(-M) all i > 

1 / T- 1 1 / t 1 



With exploitation we define an age at recruitment T r and the equations are then 

X i,t = X i-l,t-l 
X i,t = X i-l,t-l 
X o,t ' R 



The equilibrium yield Y under constant fishing mortality may be expressed as a slight 
variation of the yield per recruit expression. 

Y = R F f exp(-tM - (t - T )F) w, dt 

T r 
where w. = w (1 - exp(-kt)) 3 



The Fmax value associated with the MSY yield Ymax may then be found simply. Ymax can be 
expressed as a proportion of the unexploited recruited biomass B : 

We define B o = R / w t exp(-Mt) dt (6) 

y = Y max /B o 

and proceed to ilustrate the values of y for different values of M, K' and Tr, in Pigs. 
1-5 Superimposed on these graphs we have drawn the line y 1/2 M to illustrate the 
relationship between the supposed approximation to MSY given by equation (2) and the 
actual values calculated from the Beverton and Holt model. The details of these calcu- 
lations are given in Appendix 1. 

The essential difference between the results obtained here and those considered by 
Gulland (1971) lie in the interpretation of what is feasible, or realistic, parameter 
space. In Pigs. (1-5) where the values of y lie above the 1/2M line the age at recruit- 
ment is high and therefore the defined exploitable biomass is a small proportion of the 
total biomass. Put in another way, it is possible to obtain yields which are in excess 
of 1/2MB , as noted by Gulland, however such yields tend to occur in situations with an 
unrealistic combination of age at recruitment and natural mortality. A corollary of 
this is that they tend to occur at extremely high levels of fishing mortality and would 
in practice be probably unrealisable. 
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One alternative way of illustrating this point is by presenting the yield as a 
proportion of the total biomass of the resource Figs. 6-9- Here, as the age at 
recruitment becomes large, the yield as a proportion of the total biomass remains small 
relative to the 1 /2M line. 

Given these considerations the results are striking. For most realistic values of 
the parameters, the NSY yield as a proportion of the recruited biomass is veil below the 
1/2M level. 

The degree of overestimation varies somewhat, being greatest for the lower ages at 
recruitment and the lower values of M/K'(the ratio of mortality to growth rate). For 
some parameter combinations the yield is overestimated by up to 200?. 

There are four caveats about these results. The first is that they assume a knife 
edge recruitment. Partial recruitments will obviously effect the quantitative detail of 
the yield levels obtained. 

Secondly the equations are scaled so that the weight at the age zero is zero. Some 
stocks have "zero" weights at negative ages. Again such considerations effect the 
detail of the calulations, however it is clear that the qualitative features of substan- 
tial overestimation will remain unchanged. 

Thirdly the yield will be affected slightly if there are marked seasonal fluctua- 
tions in growth, recruitment and mortality. Preliminary studies indicate that this can 
effect the potential yield by a maximum of 20%. 

Fourthly it should be mentioned for completeness that the Beverton and Holt model, 
by assuming constant asymptotic recruitment, can be somewhat pessimistic for stocks 
which have a recruitment declining from a peak as stock size increases. Thus for stocks 
with a domed stock and recruitment relationship, if the MSY occurs at a stock level near 
the peak of the recruitment curve the results may underestimate the potential yield. 

Despite these caveats, it seems clear that the use of the 1/2 MB Q expression will 
lead, in the vast majority of cases, to an overestimation of potential yield. This is 
despite the fact that the level of yield depends on the assumption that recruitment does 
not decline with decreasing stock size. This assumption is examined in the next 
section. 

3. STOCK RECRUITMENT PROBLEMS 

It is well recognised that certain types of fish species demonstrate a decline in 
the average levels of recruitment at low stock sizes. The level of spawning stock size 
at which such declines occur is both variable, and when recruitment fluctuates, diffi- 
cult to determine. This problem is addressed in some more detail below. Here the 
possible constraints that may be chosen on escapement and how they are violated in 
certain circumstances are considered. Firstly we define the average unexploited 
spawning stock biomass as S o and the age at maturity as T m then from equation (3; 



S = R exp(-Mt) dt (7) 

T m 

and under exploitation the equilibrium spawning stock biomass S will be given by 



S = R / exp(-Mt - (t - T )F)w. dt (e) 

T m 
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if T m > T r and a slight modification if T m < T r . The ratio of the exploited to un- 



exploited stock size P is by definition 



S/S C 



In Appendix 2 the value of this ratio is tabulated for a variety of different ages 
at recruitment, mortality and growth. These results indicate that for a combination of 
high ages at recruitment and high mortality rates, harvesting at MSY will lead to a 
substantial reduction in spawning stock biomass from its unexploited level. 

Logically one should move on now to consider the necessary reduction in yield that 
would ensure that the SSB was not substantially reduced. However, because the Beverton 
and Holt model has such a flat yield curve, the level of reduction required is, in many 
cases, rather slight. By contrast the necessary reduction in fishing mortality is 
substantial. A typical case is illustrated in Fig. 10. 

For the situation illustrated in Fig 10a, the MSY is around 11# of the unexploited 
recruited biomass, a reduction in the fishing mortality from its Fmax valae of around 
0.52 to about half that level will approximately double the equilibrium SSB, but will 
only reduce the catch by about 1 0%. However it should be emphasised that this is a 
rather artificial result. In Fig 10b the relation between the equilibrium spawning 
stock and catch level is illustrated. It indicates that very small fluctuations in 
catch around the MSY level will have large effects on the spawning stock. Accordingly 
as will be shown below, where recruitment varies randomly these deterministic results 
are effectively meaningless. 

In Table 1 we present a summary of the results obtained so far, for typical values 
of mortality, growth rate, and ages at maturity and recruitment. The implications of 



Table 1 



Age at recruitment 



1 



Mortality Rate 
.2 .4 .6 



Mortality Rate 
.2 .4 .6 



Mortality Rate 
.2 .4 .6 



Y/Brec .042 .067 .090 
Y/BTot .042 .067 .090 
S/So .293 -307 .274 



.049 .085 .125 
.049 .085 .124 
.267 .261 .214 



.058 .112 .180 
.057 .107 .162 
.245 .219 .179 



Y/Brec .055 .084 .110 

0.4 Y/BTot .055 .084 .110 

.269 .269 .252 



.068 .116 .167 
.067 .115 .163 
.227 .206 .162 



.085 .163 .256 
.082 .148 .213 



S/So .269 .269 .252 .227 .206 .162 .196 .162 .113 

MSY levels expressed as a proportion of unexploited recruited biomass (Y/Brec)and as a 
proportion of the total biomass. (Y/BTot) The levels of equilibrium spawning biomass 
under harvesting.(S/So) are also shown. The age at maturity for these results is 2 
years . 



these results may briefly be summarised. In most of the feasible biological situations 
the maximum yield that can be taken from a stock is well below the level given by 
1/2MB Furthermore and in particular for the high mortality rates, there is a further 
need to reduce the estimates of potential yield to ensure that the spawning stock 
biomass is not reduced to a level where recruitment may be expected to decline. 

These results are also subject to a caveat. In the analysis, we have made the 
optimistic assumption that the life history occurs in a sequence in which recruitment to 
the adult stock is followed by spawning and then harvesting. This means when the age at 
recruitment is greater than or equal to the age at maturity, that the spawning stock 
level is kept at a maximum. If harvesting occurs prior, to or during spawning, the 
effect it has on SSB will be somewhat greater. This is particularly true for stocks 
with the higher mortality rates. 

So far the results have been obtained on the assumption that recruitment is constant. 
Ve now proceed to examine the problems raised by variable recruitment. 

4. THE PROBLEM OF RECRUITMENT VARIATION 

A number of recent studies have examined the variation in recruitment for a variety 
of fish stocks. The basic work is that of Hennemuth, Brown and Palmer (1980), and 
subsequent analysis has little extended their conclusions. These conclusions are, that 
recruitment is highly variable for some stocks and tends to possess a skewed frequency 
distribution, the log normal being a reasonable empirical fit. 

In Table 2 we present the results of an analysis of a variety of data sets in which 
the recruitment variation is presented as the variance of the log normal distribution. 
This statistic has the advantage of being comparable across data sets with different 
absolute levels of recruitment. In the subsequent analysis we use the square root of this 
statistic, that is, the standard deviation of the log-recruitment, which we denote by 0** 

The data sets used are, in a sense, self selecting. This is because to obtain 
information on recruitment variation, a knowledge of the catch at age or length and the 
natural mortality rate is required. Hence the data are from well studied stocks, for 
which cohort or virtual population analyses have been performed. To this extent they 
cannot be considered typical. Furthermore the values of the standard deviation given 
are only a very approximate guide to the level of variability: recruitment estimates for 
at least 20 successive years would be necessary to estimate the S.D. to within a factor 
of two, even if the annual recruitments were mutually independent. In practice years of 
good or poor recruitment tend to occur together, and so the margin of error in estim- 
ating the standard deviation is even greater than this. In all that follows we shall 
assume that successive annual recruitments are mutually independent provided that the 
spawning stock biomass is not reduced below some specified minimum level. This assump- 
tion will tend to result in an underestimation of the effects of variability. 

5. THE IMPLICATIONS OP RECRUITMENT VARIATION FOR ESTIMATES OF POTENTIAL YIELD 

There are two basic problems posed by recruitment variation, in estimating the 
potential yield of fish stocks. The first involves the problem of declining recruitment 
with stock size. In this situation it is possible that stochastic variation in recruit- 
ment will, with high probability, ensure that the SSB will be reduced below some fixed 
value. This can occur even though in the deterministic case, the equilibrium SSB is 
above that value. 

The second problem is that in the process of survey and assessment of potential 
yield, the survey is effectively sampling from an underlying probability distribution of 
the biomass of the fish stock. Hence if the biomass is at an uncharacteristically high 
level, potential yield could be overestimated, or conversely, underestimated, if the 
biomass is uncharacteriscally low. 
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In analysing these tiro problem* we have related them, by investigating, using Monte 
Carlo techniques, the whole process of survey, assessment and harvesting. In this way we 
can assess the probability that the stock will be reduced to some constraint level under 
different levels of harvest. 

The basic question that can be answered from the Monte Carlo investigations is Vhat 
is the probability that the spawning stock biomass will drop below some specified 
escapement level in a fixed period of time? 1 However to pose the question in the 
context of assessing potential yields we have asked the question. 'What modification of 
the basic catch level will ensure that this probability is less than some specified 
amount? ' 

In order to gain some insight into the process, prior to considering the stochastic 
system, it is useful to consider the dynamics of the deterministic case. Once 
harvesting begins, the stock biomass will follow a pattern of the type illustrated in 
Figure 11. Harvesting at a constant fishing mortality will give catches which are 
initially large, but which decline rapidly as the population moves towards the new 
equilibrium. By contrast, the operation of a constant catch quota will require 
initially a small fishing mortality, which will necessarily be increased towards a 
maximum as the equilibrium biomass is approached. The approach to equilibrium in this 
case will be slow relative to the constant fishing mortality strategy. 

The biological characteristics of the resource also have an effect on the speed of 
approach to the equilibrium. The mortality rate, and to a lesser extent the von 
Bertalanffy growth coefficient, determine the time scale of the response to exploit- 
ation. When both are high the approach to equilibrium is rapid, when both are low, it 
may take two or more decades. 

Superimposed on this deterministic process are the two stochastic factors caused by 
recruitment variation. One affects the initial biomass and hence, in absolute terms, 
the level of harvest attempted; the other directly affects the population dynamics as 
exploitation progresses. We now explore the consequences of this variation. 

In a system where recruitment is stochastic, instead of an unexploited equilibrium 
stock size, there is a probability distribution of recruited population biomass. The 
expected value of that distribution is 



E (B ) = R r exp(-Mt) w t dt (9) 

T r 

where R is the mean recruitment rate, and its variance is: 



var (B Q ) = v / exp(-2Mt) w dt (10) 

where v is the variance of the recruitment rate. We ignore here possible year-to-year 
variations in the growth and mortality rates, in practice these variations exist, but 
have much less effect these recruitment variations. 

In a similar way the expected values of the spawning stock biomass can be written 
down once the age at sexual maturity is specified. It is then possible to specify some 
escapement constraint* In all the results that follow, an escapement level of 20% of 
the expected unexploited spawning stock biomass is used. This is not a conservative 
figure, but it represents a lower limit where recruitment declines might be expected to 
be observable. The time scale that has been used is similarly arbitrary, although its 
choice has been guided by common sense. We have chosen a twenty year period in which to 
investigate the probability that the escapement will fall below the 20% level. 
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In presenting the results of this analysis, we have calculated the appropriate level 
of catch, .that irill ensure that the probability that the SSB falls below 20% of its 
unexploited level is less than 0.1* 

As we have already indicated. The different strategies of a constant catch and a 
constant fishing effort have rather different dynamic consequences and hence it is 
necessary to present the results in two different ways. 

5-1 Constant Catch 

The calculations here are presented as modifications of the basic deterministic 
results presented earlier. In Figs (12 -15) the results of the Monte Carlo investi- 
gations have been analysed so that the catch level that can be taken as a proportion of 
the initial surveyed biomass is presented as a function of mortality, growth and the age 
at recruitment. This catch level is the maximum one that ensures that the probability 
that the SSB will be reduced to less than 20% of its expected level in 20 years is less 
than 0.1. 

In all the calculations we have used an age at sexual maturity of two. This will be 
either accurate, or optimistic, for most fish stocks. However for fish stocks with 
mortality rates in excess of 0.?, the results will be pessimistic as such stocks 
typically mature at an age less than two. 

The results indicate that for higher levels of recruitment variation it may be 
necessary to further adjust the catch level downwards to ensure that SSB remains above 
its constraint level. Table 3 summarises these results. 



5*2 Constant Fishing Mortality 

The results obtained for a strategy involving constant fishing mortality are not 
directly comparable to those that can be presented for constant catch. Instead the 
process needs to be considered in two stages. Figs. 16-19 illustrates the reductions 
in fishing mortality from Fmax that are required to satisfy the constraint that the 
probability that the SSB will fall below this level within 20 years should be less than 
0.1. 

The results are illustrated for different values of recruitment variability, 
Mortality, growth and age at recruitment. Associated with each set of these parameters 
and constraints is a fishing mortality rate. If applied to the stock this fishing 
mortality rate produces a catch level that reduces as the equilibrium is approached. 

Figs. 20-23 illustrate these changing expected catch levels. The results show that 
the constraints on the SSB involve substantial sacrifices in the initial catch compared 
with those obtainable from fishing at Fmax, but no or only modest sacrifices in the 
longer term catch. 

In Figs. 24 - 25 the approximate confidence regions around some of these expect- 
ations are presented. They indicate that where the recruitment variation is large, 
variation in catch is so great that the practicality of operating a constant fishing 
Mortality is questionable. 

In Table 4 the results of the calculations on the strategy of constant fishing 
Mortality are summarised. 

This expected level is defined in equation (9) and may best be considered as the 
long term average abundance of the fish stock in its unexploited state. 
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Table 3 
Constant catch strategy 

Age at recruitment / sexual maturity 

1 2 

Mortality Rate Mortality Rate 

.2 .4 .6 .2 .4 .6 

K' 

unconstrained .049 .085 .125 .058 .112 .180 

0.2 cr-0.4 .049 .082 .104 .058 .104 .149 

<r-1.0 .049 .077 .084 .058 .092 .104 

unconstrained .068 .116 .167 .085 .163 .256 

0.4 (T-0.4 .068 .101 .129 .085 .139 .200 

(7-1.0 .068 .086 .092 .085 .106 .110 

MSY levels expressed as a proportion of the mean unexploited recruited biomass subject 
to the constraint that the probability spawning stock biomass falling below 20% of its 
initial level within 20 years after the start of fishing at a constant catch level is 
less than 0.1 , with two levels of recruitment variability, cr (see text). The uncon- 
strained MSY is shown for comparison. 
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Table 4 



Constant fishing mortality 



Age at Recruitment 



Mortality Rate 



Mortality Rate 



K 
0.2 



unconstrained 
<r - 0.4 
<? - 1.0 



.2 .4 .6 

.042 .067 .090 
.042 .067 .090 
.041 .061 .072 



.2 .4 .6 

.058 .112 .180 

.058 .108 .173 

.055 .092 .121 



0.4 



unconstrained 

<T - 0.4 



.055 .084 .110 
.055 .084 .109 
.052 .069 .075 



.085 .163 .256 
.084 .163 .236 
.073 .112 .131 



MSY under a constant 7 strategy: the figures represent the maximum long-term yields (as 
a proportion of average unexploited biomass) subject to the constraint that the proba- 
bility of the spawning stock biomass being reduced to below 20% of its average 
unexploited level within 20 years is less than 0.1. Two values of recruitment varia- 
bility, a , are used and the unconstrained MSYs shown for comparison. Age at sexual 
maturity 2 years. 



6. THE IMPLICATIONS FOR-SURVEY DESIGN 

It is immediately obvious that there is one major implication of the analysis. 
Clearly future surveys should 9 in addition to estimating biomass, attempt at least to 
obtain sufficient samples to estimate the parameters of growth and mortality. 
Information on recruitment variability is also important, and where a species is rela- 
tively long lived, some idea of this can be obtained from the same age/length samples 
required to estimate the other parameters. 
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7- THE CALCULATION OF POTENTIAL YIELD FROM A SURVEY 

The whole of the foregoing analysis can now be encapsulated in a series of calcu- 
lations and checks that are needed to estimate a reasonable range for the potential 
yield of a fish resource. 

(1) Given information on the biomass, mortality and growth, a range of estimates of the 
MSY can be obtained for a chosen age of recruitment. This can be done using Pigs. 1 
-5 or Appendix 2 Table 1. 

(2) A check should then be made to see if such a level of catch would be likely to 
reduce the SSB below some target level. This can be done using Appendix 2 Table J. 

(j) The level of recruitment variation to be expected from nch a species should then be 
considered. Often this will need to be done by analogy with similar species in 
other f hopefully, similar areas. 

(4) Given a level or range of recruitment variation the modification in catch level fog 
a constant catch strategy can be calculated. This can be done using Figs. 12-15 

A parallel calculation can be made for a constant fishing mortality strategy. The 
short term and long term catch levels to be expected can then be obtained by 
inspecting Figs. 20 - 2?. 

(5) The potential yield of the stock can then be presented as a range calculated for the 
different strategies of constant catch and constant fishing mortality. Where con- 
siderations of a decline in spawning stock size are irrelevant these two strategies 
will lead to similar results. 

It should be emphasised that the estimates obtained in this way are approximate. 
Indeed even where parameters can be estimated directly, if the full range of uncertainty 
is considered the yield will usually estimated only to within a factor of two or more. 



These figures have been constructed, as already discussed, on the basis that the 
probability of a SSB decline to below the constraint level, is less than 0.1, in a 
20 year period. 
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F1g. 2 
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Fig. 8 
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F1g. 9 
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Fig. 11 
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Fig. 12 

Potential yield as a percentage of BO subject to the constraint 
that the spawning stock biomass remains above 20% of the unex- 
ploited level in a 20 year period with a probability of 90%. 
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F1g. 13 

Potential yield as a percentage of BQ subject to the constraint 
that the spawning stock biomass remains above 2O% of the unex- 
ploited level in a 2O year period with a probability of 9O%. 
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Fig. 14 

Potential yield as & percentage of BQ subject to the constraint 
that the spawning stock biomass remains above 2O% of the unex- 
ploited level in a 2O year period with a probability of 9O%. 
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F1g. 15 

Potential yield as a percentage of B o subject to the constraint 
that the spawning stock biomass remains above 20% of the unex- 
ploited level in a 20 year period with a probability of 90%. 
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Flg. 16 

Reductions in fishing mortality from Fmax required to satisfy 
the constraint that, with stochastic recruitment, the probabr- 
ility that the SSB falls to below 20% of its average initial 
level is less than 0.1, where cr * 0.4 and 1.0 as indicated. 
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Fig. 17 

Reductions in fishing mortality from Fmax required to satisfy 
the constraint that, with stochastic recruitment, the probabr 
ility that the SSB falls to below 20% of its average initial 
level is less than 0.1, where cr - 0.4 and 1.0 as indicated. 
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Fig. 18 

Reductions in fishing mortality from Fmax required to satisfy 
th<2 constraint that, with stochastic recruitment, the probabi- 
lity that the SSB falls to belovr 2O% of its average initial, 
level is less than 0.1, where cr O.4 and 1.0 as indicated. 
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Fig. 19 

Reductions in fishing mortality from 'Fmax required to satisfy 
the constraint that, with stochastic recruitment, the probab- 
ility that the SSB falls to belovr 2O% of its average initial 
level is less than O.1, where a * 0.4 and 1.O as indicated. 
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APPENDIX 1 

Derivation of maximum yield calculations 

The derivation here is essentially the same as that given by 
Beverton and Holt (1957). 

For a recruitment rate, R, at age zero, the number of fish in 
the unexploited state aged t is Re""*^ where M is the natural 
mortality rate. We assume that the growth follows a von Bertalanffy 
formula, with a zero weight at age zero, ie: 

w t = w^ (1 - e" kt ) 3 (1) 

where w t is the weight at age t and w^ the asymptotic average 
weight. This implies that the biomass flux at age t is: 

Rv^e-^d - e' kt ) 3 (2) 

and so the total biomass is: 

Rw~ /" e" tM (1 - e~ kt ) 3 dt (3) 

which is equal to: 

6K 3 Rw /{M(M + K) (M + 2K) (M + 3K) } (4) 

OO 

If we define an age at recruitment, t r/ then the recruited 
biomass is: 

Rw^ /" e'^d - e" Kt ) 3 dt (5) 

which is equal to: 

Rw e -V ( 1 3e ~ ktr + 3e ~ 2kt;r - e " 3ktr } (6) 
Rw co e r { M M + K + M + 2K M - 3K } (6) 

The equlibrium recruited biomass under constant fishing mortality 
F is thus: 

- RW e-V { i - 3e" kt r . 3e" 2kt r e" 3kt r } (7) 
Z Z+K Z+2K~Z+3K 

where Z = M + F. 
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The age at recruitment t that should be used n this formula is 
that referred to the theoretical age at zero weight extrapolated 
from the von Bertalanffy growth curve. If this is non-zero , then 
the value for t to be put in .this formula will not be the actual 
age at recruitment. 

If the fishing mortality rate is constant, then the catch rate 
will be: 

FRw e'V { 1 - 3e ~ ktr + 3e" 2kt r _ e~ 3kt r } (8) 
* KW e ri Z Z + K Z + 2K Z + 3K ' 

If the catches occur discretely , that is, are confined into 
short annual seasons, during which the growth and natural mortality 
rates can be assumed to be negligible in comparison, then the 
annual catch is: 



- e ~ (z+k) 



e -3kt 

6 r } (9) 



1 _ e -(z+3k) 

In compiling figures 1-9 and the tables in Appendix 2, formula 
(8) rather than (9) has been used. 

Equation (8) can also be expressed in terms of the length at 
recruitment. If c is the length at recruitment as a proportion 
of the average maximum length, then the yield is: 

FRw (1 - c) M/K { 1 - 3(1 ~ c > + 3(1 - c) 2 _ d - c) 3 } (10) 
J?KW ~< I C) 1 Z Z + KZ + 2K Z + 3K' 

To find the maximum yield for a given length or age at recruitment, 
we simply maximise (8) or (1O) over F. If c is less than the 
eumetric length c 1 (the length at which growth and natural mortality 
exactly balance) given by: 



then the yield is obtained for a finite value of F. There is 
no simple formula for F , but it can easily be evaluated numerically 
or graphically by plotting a yield versus F curve, such as the one in 
Figure 1Oa. Such curves are often rather flat-topped, with a whole 

range of F values either side of F mav all giving much the same yield. 

max 

If c is greater than c 1 , then the yield increases asymptotically 
to a maximum value as F increases. 



The theoretical maximum yield obtainable in terms of the recruited 
unexploited biomass, B O , from these formulae is MB Q , but this will 
only be achieved at very high ages or lengths at recruitment and in 
practice the yields will be much less. Figures 1-5 show the maximum 
yields as a proportion of B Q for various parameter values. 

The maximum yield as a proportion of total initial biomass is given 
in figures 6-9 , for various parameter values. The maximum possible 
value is achieved by eumetric fishing, ie. where c = c 1 and F is 

very large. The yield is then a MB , where: 

max o 



a mav 9 U - C ' ) "( 1 - C 1 + "S^ ' * ) U2) 

1UCLA, ^T" 17 

where y = M/K. 

If c > c 1 , the yield is aMB Q , where: 

a - c 3 (1 - c) u (1 + y) (2 + y) (3 + y) (13) 
6 

If c < c ' , there is no simple exact formula for a, but it can 
be approximated to within +/- 1O% by the formula: 

a a max (1 + c/c , } (14) 

In most cases c will be less than c 1 , often much less. 

The spawning stock biomass corresponding to a given level of 
fishing mortality is: 

where t is the age at sexual maturity. If t > t , which will 
m m r 

nearly always be the case in practice r then the spawning stock 
biomass is equal to: 



Variable recruitment 

If we suppose that the recruitment rate is a random variable, with 
successive annual recruitments independently and identically 
distributed about a mean R with a coefficient of variation s 
(ie a variance of R 2 s 2 ) then from equation (5) the variance of 
the recruited biomass is: 
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R 2 s 2 w / e~ 2tM (1 - e~ Kt ) 6 dt (17) 

The integral can be expanded in an analogous form to (6) . Likewise 
analogous formulae to (6) ~ (10) can be obtained for the variances 
of the different quantities* 

A statistic of interest is the coefficient of variation of the 
catch and biomass under fishing at constant F r which can easily 
be calculated from these formulae. In practice, the parameter s, 
will not be known with any precision, and an approximate version 
of the formula for the co-efficient of variation will suffice. 
The approximation: 

c.o.v. = - s /Z (18) 

holds quite well for all parameter values and can be applied to 
the catch, the recruited biomass, or the spawning stock biomass. 

The variation in recruitment is expressed in Table 2 in terms of 
a slightly different statistic, the standard deviation of the 
log-recruitment, a. If the log-recruitment is normally distributed 
then s is related to or through the formula: 



(19) 

Equation (18) should be used cautiously for three reasons. Firstly, 
recruitment in fish stocks tends to have a very skewed distribution 
and so the confidence limits on the catch and biomass will be 
asymmetrical. Secondly, in most fish stocks, successive annual 
recruitments are not mutually independent, but show a strong 
serial correlation. Even relatively weak serial correlation will 
invalidate (18). Thirdly, the sample sizes will generally be too 
small to estimate s or a to any degree of precision, and this will 
be exacerbated by any serial correlation that may be present. 
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APPENDIX 2 



Table 1* Maximum sustainable yields expressed as a proportion of the 

unexploited recruited biomass, for different ages of recruitment, 
natural mortality rates, and values of the von Bertalanffy growth 
parameter, K. 



AGE AT RECRUITMENT t 



AGE AT RECRUITMENT! 3 





N 


















.1 


.2 


.3 


.4 


.3 


.6 


.7 


.8 


K 


















.1 


.021 


.033 


.043 


.056 


.067 


.078 


.088 


.099 


.2 


.027 


.042 


.033 


.067 


.078 


.090 


.101 


.112 


.3 


.032 


.049 


.063 


.076 


.088 


.100 


.112 


.123 


.4 


.036 


.033 


.070 


.084 


.097 


.110 


.'22 


.134 


.5 


.03? 


.060 


.077 


.091 


.105 


.118 


.131 


.143 


AGE AT 


RECRUITMENT! 


1 














M 


















.1 


.2 


.3 


.4 


.5 


.6 


,7 


.8 


K 


















.1 


.023 


.037 


.052 


.068 


.083 


.103 


.122 


.143 


.2 


.030 


.049 


.067 


.085 


.103 


.125 


.147 


.171 


.3 


.036 


.039 


.080 


.101 


.123 


.146 


.171 


.197 


.4 


.041 


.068 


.092 


.116 


.141 


.167 


.195 


.224 


.5 


,046 


.076 


.103 


.130 


.158 


,18- 


.218 


.230 


AGE AT 


RECRUITMENT: 


2 














M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.024 


.043 


.063 


.085 


.111 


.141 


.173 


.213 


.2 


.034 


.058 


.084 


.112 


.144 


.180 


.220 


.265 


.3 


.042 


.072 


.104 


.138 


.176 


.218 


.265 


.316 


.4 


.048 


.085 


.123 


.163 


.207 


.256 


.309 


.368 


.3 


.054 


.097 


.141 


.187 


.238 


.293 


.333 


.418 





M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.027 


.049 


.073 


.107 


.143 


.190 


.243 


.304 


.2 


.038 


.069 


.104 


.145 


.193 


.249 


.313 


.383 


.3 


.048 


.088 


.132 


.183 


.241 


.307 


.381 


.438 


.4 


.036 


.103 


.159 


.219 


.287 


.362 


.441 


.522 


.5 


.064 


.121 


.183 


.233 


.329 


.409 


.492 


"".577 


A6E AT 


RECRUITMENT} 


4 














N 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.029 


.036 


.090 


.132 


.189 


.248 


.322 


.39? 


.2 


.043 


.081 


.128 


.184 


.251 


.326 


.406 


.488 


.3 


.034 


.103 


.164 


.234 


.311 


.392 


.476 


.562 


.4 


.064 


.127 


.197 


.276 


.359 


.444 


.532 


.622 


.3 


.073 


.143 


.226 


.310 


.397 


.486 


.376 


.668 


A6E AT 


RECRUITMENT: 


5 














M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.032 


.064 


.106 


.161 


.230 


.306 


.387 


.472 


.2 


.048 


.095 


.154 


.226 


.305 


.388 


.474 


.562 


.3 


.061 


.123 


.197 


.278 


.363 


.451 


.341 


.632 


.4 


.072 


.148 


.231 


.317 


.407 


.498 


.390 


.684 


.3 


.081 


.166 


.254 


.345 


.438 


.531 


.626 


.722 



Table 2 : As Table 1 , but with yields expressed as a proportion of total 



AGE AT RECRUITMENT: 


unexploited 




biomass. 




N 


















.1 


.2 


.3 


.4 


.3 


.6 


.7 


.8 


K 


















.1 


.021 


.033 


.043 


.036 


.067 


.078 


.068 


.099 


.2 


.027 


.042 


.055 


.067 


.078 


.090 


.101 


.112 


.3 


.032 


.049 


.063 


.076 


.088 


.100 


.112 


.123 


.4 


.036 


.055 


.070 


.084 


.097 


.110 


.122 


.134 


.5 


.039 


.060 


.077 


.091 


.105 


.118 


.131 


.143 



AGE AT RECRUITMENT: 1 



.2 



.3 



.4 



.5 



.7 



.8 



K 


















.1 


.023 


.037 


.052 


.068 


.083 


.102 


.121 


.141 


.2 


.030 


.049 


.067 


.085 


.104 


.124 


.145 


.167 


.3 


.036 


.059 


.080 


.100 


.122 


.144 


.167 


.191 


.4 


.041 


.067 


.091 


.113 


.139 


.163 


.188 


.214 


.5 


.045 


.075 


.102 


.128 


.155 


.181 


.208 


.236 


AGE AT 


RECRUITMENT: 


2 














H 


















.1 


.2 


.3 


.4 


.3 


.6 


.7 


.8 


K 


















.1 


.024 


.042 


.062 


.083 


.107 


.132 


.159 


.187 


.2 


.034 


.057 


.081 


.107 


.134 


.162 


.191 


.220 


.3 


.041 


.071 


.099 


.129 


.159 


.189 


.219 


.249 


.4 


.048 


.082 


.115 


.148 


.181 


.213 


.244 


.273 


.5 


.033 


.093 


.130 


.165 


.200 


.233 


.264 


.293 



A8E AT RECRUITMENTS 3 





M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.027 


.048 


.072 


.098 


.127 


.135 


.183 


.209 


.2 


.037 


.066 


.096 


.126 


.156 


.185 


.211 


.233 


.3 


.046 


.082 


.116 


.149 


.180 


.207 


.230 


.243 


.4 


.054 


.095 


.132 


.167 


.197 


.222 


.239 


.248 


.5 


.060 


.105 


.145 


.180" 


.209 


.22V 


.241" 


.244 


A6E 


AT RECRUITMENT: 


4 














M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 




K 


















.1 


.029 


.034 


.081 


.110 


.138 


.164 


.183 


.193 


.2 


.041 


.074 


.106 


.137 


.163 


.182 


.191 


.191 


.3 


.051 


.090 


.125 


.155 


.176 


.186 


.187 


.180 


.4 


.059 


.102 


.138 


.164 


.178 


.182 


.177 


.166 


.5 


.064 


.1)1 


.146 


.167 


.175 


.173 


.165 


.131 


AGE 


AT RECRUITMENTi 


5 














M 


















.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


K 


















.1 


.031 


.039 


.088 


.116 


.140 


.133 


.135 


.147 


.2 


.045 


.080 


.112 


.137 


.150 


.131 


.143 


.130 


.3 


.055 


.095 


.126 


.143 


.146 


.140 


.127 


.111 


.4 


.062 


.104 


.130 


.140 


.137 


.127 


.112 


.093 


.5 


.066 


.108 


.129 


.134 


.127 


.114 


.098 


.081 
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